Tuning the charge-transfer energy in hole-doped cuprates 
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Chemical substitution, combined with strain, allows the charge-transfer energy in hole-doped 
cuprates to be broadly tuned. We theoretically characterize the structural and electronic properties 
of the family of compounds .R2CUO2S2, constructed by sulfur replacement of the apical oxygens and 
rare earth substitutions in the parent cuprate La2Cu04. Additionally, the enthalpies of formation 
for possible synthesis pathways are determined. 



In weakly correlated systems, analytic and computa- 
tional tools, such as density functional theory (DFT), 
have achieved a level of control where one can reason- 
ably contemplate materials design. For recent work in 
this upcoming field, see Refs. [1-5]. The corresponding 
methods for strongly-correlated systems are rapidly ad- 
vancing [6-11]. Materials design of correlated systems 
is a valuable route to test our physical understanding, 
and should the materials be synthesized, provides exper- 
imental feedback to accelerate development of theoretical 
methods. 

The high-temperature cuprate superconductors are 
classic strongly-correlated systems which exhibit the 
highest transition temperatures known. Despite intense 
theoretical and experimental study, we have limited first- 
principles understanding of the chemical parameters con- 
trolling these transition temperatures. Based on an ab 
initio study of the cuprates [12], we proposed that the 
charge-transfer energy — e p controls the strength of 
correlations [13] and therefore tunes the maximum super- 
conducting transition temperature T Cjmax . Starting with 
the most correlated cuprate La2Cu04, reducing — e, p 
enhances T Cimax . In contrast, Sakakibara, et. ai, pro- 
posed that reducing the admixture of the Cu 3d z 2 or- 
bital with the in-plane 3d x 2_ y 2 orbital enhances T Cjmax , 
a mechanism dubbed "orbital distillation" [14,15]. These 
studies pose a challenge for materials design to test the 
theoretical proposals. 

In this work, we design a family of copper oxysulfides to 
tune the charge-transfer energy of La2Cu04 (LCO) and 
test the proposed theories. We show that substitution of 
the apical oxygens by sulfur produces sufficiently large 
effects on td — £ p and orbital distillation to warrant in- 
vestigation. The compounds are valuable for differentiat- 
ing the two theories since, barring competing orders, the 
proposed mechanisms predict opposite shifts in T C]lnax in 
the oxysulfides relative to the LCO baseline. However, 
checks on structural stability and the thermodynamics 
of synthesis pathways indicate epitaxial stabilization is 
likely necessary for realization of the oxysulfides. 

We choose the T-type layered perovskite La2Cu04 
(LCO) as our structural starting point. This single- 
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FIG. 1. Scatter plot of the charge-transfer energy (ed — e p ) and 
the apical atom distance from the CUO2 plane (d ap j ca i) for the 
family of hypothetical copper oxysulfides R2CUO2S2 (squares) 
and the actual compound La2Cu04 (circle). The oxysulfides 
are formed by the substitution of apical oxygens by sulfur 
(arrows in inset figure of crystal structure), and the structure 
obtained by full relaxation in DFT. Chemical substitution 
allows us to span nearly the entire range of charge-transfer 
energies (shaded bar) found in the cuprates. 

layered cuprate is well-studied, supports a broad range 
of doping, and exhibits the largest charge-transfer energy 
of all the hole-doped cuprates. Guided by the intuition 
that larger, more covalent, apical ions would decrease 
£d — e pi we replace the apical oxygens in LCO with sulfur 
to form the copper oxysulfide La2Cu02S2- Due to the 
larger ionic radii of sulfur relative to oxygen, we expect 
the LaS charge reservoir layers to be crowded. To com- 
pensate for this large sulfur in the charge reservoir layers, 
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FIG. 2. Fully relaxed conventional cell lattice parameters for 
R2CUO2S2 family of compounds (squares) as compared to 
La2Cu04 (dot). Sulfur substitution increases both the a and 
c lattice parameters. The decrease in the cation's ionic radii 
from (La— >Y— »Lu— >Sc— >Gsl) accounts for the contraction of 
a and c within the family of oxysulfides. The dotted line is a 
guide to the eye. 



we explore the effect substitution of the large La ion with 
smaller trivalent ions R =Y. Lu, Sc, Ga. We structurally 
relax, within a single unit cell, the entire family of cop- 
per oxysulfides i?2Cu0 2 S 2 , using the Vinet equation of 
state [16]. We used the VASP [17-19] density functional 
software package with PAW potentials [20, 21]. Finally, 
we extract td — e p by performing a full-potential calcu- 
lation on the relaxed structure using Wien2K [22], then 
downfolding [23] to orbitals constructed in the manner 
described in Ref. [24]. 

The charge-transfer energies for the parent compound 
La2Cu04 and the family of hypothetical oxysulfides 
i?2CuC>2S2 are shown in Fig. 1. The substitution of 
sulfur alone, leaving the La ion unmodified, decreases 
the charge-transfer energy by 0.6 eV, a significant reduc- 
tion given that 64 — e p spans a range of 1.4 eV across all 
cuprate families. If no competing electronic or structural 
orders are introduced, we predict La2Cu02S 2 to have 
a higher superconducting transition temperature than 
LCO (which has a T c max of 38 K). Further substitution 
of La by trivalent cations of varying ionic radii allows us 
to span nearly the entire range of charge-transfer ener- 
gies found in the cuprates. We find that the apical atom 
distance d a picai is increased by as much as 0.2 Ain the 
oxysulfides, due to the large size of sulfur. 

The structural effects of the sulfur and rare-earth sub- 
stitutions can be rationalized entirely due to trends in 
ionic radii, as shown via a scatter plot of the a- and c-axis 
lattice constants of the proposed compounds (Fig. 2). 
The ionic radius of sulfur is greater than that of oxygen 
by 0.4 A, a significant amount. Apical sulfur substitution 
alone enhances the c-axis length by over 1.5 Arelative 
to the parent compound LCO, in agreement with the 
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FIG. 3. The effect of uniaxial c-axis strain on the charge- 
transfer energy (circles) for the actual compound La2Cu04 
(top) and the oxysulfide La2Cu02S2- The dotted line shows 
the equilibrium d ap icai- In LCO, uniaxial compression in- 
creases ed — e p and consequently suppresses T c , as observed 
experimentally [25], while we predict pressure should hardly 
affect T c in the oxysulfide. To disentangle crystal fields from 
the Madelung potential, we show the effect of tweaking d ap i ca i 
with all other atoms fixed (squares). The effects due to the 
Madelung potential, roughly quantified by the difference be- 
tween the two curves, is the same order of magnitude as crys- 
tal fields. 



fact that there are four sulfur layers in the conventional 
cell, giving Ac ps 4 x 0.4 A. This chemically induced 
strain may be beneficial for superconductivity since it 
is known uniaxial compression strongly suppresses T c in 
La2Cu04 [25], a point we will subsequently address. In 
contrast, the in-plane constant exhibits only a slight ex- 
pansion by 0.2 A, due to the restraining effect of the stiff 
Cu-0 bonds. The subsequent rare-earth substitutions in 
the oxysulfides drive a clear structural trend, with both 
lattice parameters contracting in parallel with the de- 
creasing ionic radius of the rare-earth ion R. The same 
trend is experimentally observed in the T'-type -R2CUO4 
family [26]. 

In order to disentangle the structural contributions 
driving variations in the charge-transfer energy, we track 
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FIG. 4. Octahedral rotations in SC2CUO2S2 as shown by (a) 
a view along the a-axis, and (b) a section of the Cu02 plane. 
We structurally relaxed a 2 x 2 x 1 supercell with no symmetry 
constraints. While there is no out-of-plane buckling of the 
Cu-0 bonds, the small Sc ion causes strong in-plane rotations 
of the octahedral cages (denoted a°a°Cp in Glazer notation). 
The compound remains in the tetragonal T-type phase. 

the change in — e p as we (1) vary the c-axis parame- 
ter, modeling uniform uniaxial strain and (2) vary d a picai 
alone with all other atomic positions fixed. The differ- 
ence in the resultant — e p between these two structural 
configurations provides a rough proxy for the magnitude 
of effects due to the Madelung potential as compared 
to the local crystal fields generated in part by the api- 
cal atom. For the parent compound LCO (Fig. 3 top), 
c-axis compression enhances e^ — e p . This increase in cor- 
relation strength provides a natural explanation for the 
observed suppression of superconductivity with uniaxial 
pressure [25] . Modifying d a picai alone has a similar, albeit 
smaller, effect. The difference between the two curves is 
of the same magnitude as the total shifts in ed — e P , so 
we conclude that crystal fields and Madelung potentials 
have similar energy scales, making first-principles calcu- 
lations necessary for quantitative insight. In comparison, 
strain has the opposite effect in the hypothetical com- 
pound La2Cu02S2 (Fig. 3 bottom), with c-axis compres- 
sion driving a reduction in e^ — e p . Thus, T cmax should be 
observed to trend in the opposite direction as compared 
to the parent compound. 

Finally, we check the structural and thermodynamic 
stability of the hypothetical oxysulfides. Determining 
the lowest-energy structure of a general multi-element 
compound is extremely challenging, and the checks we 
perform are by no means exhaustive. However, guided 
by knowledge of typical competing phases and structural 
distortions, we show the .R2CUO2S2 family passes several 
basic criteria for stability. 

In La2CuC>4, when the large La ion is replaced by the 
smaller Nd and Pr atoms, the oxygens in the rock-salt 
LaO layers in the T-type structure tend to be unstable 
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TABLE I. Synthesis pathways for various cuprate oxysulfides 
based on substitution of sulfur for both (top block) or only 
one (middle block) of the apical oxygens in i?aCu04. Ener- 
gies in kj/mol and volumes in kJ/mol/GPa. Since the en- 
ergies of formation (AE = £ pro ducts - .Eteactants) are pos- 
itive, none of these pathways appear favorable at ambient 
conditions. However, high-pressure synthesis will help sta- 
bilize these pathways, since the majority of volume differen- 
tials (AV = ^products — Keactants) are negative. We bench- 
mark our method against the standard synthesis pathway 
for LaaCu04, shown on the last line. Surprisingly, AE is 
+28 kj/mol, so either DFT systemmatically overestimates 
enthalpies (which means the actual enthalpies for our hypo- 
thetical compounds are smaller, in our favor), or we must add 
a bi-directional uncertainty of ±30 kj/mol to the computed 
enthalpies. Additionally, positional entropy of the apical S in 
the half-substituted .R2CUO3S compounds should also assist 
in synthesis. 



towards rearrangement into the fluorite structure, form- 
ing the T'-type structure associated with the electron- 
doped cuprates [27]. We check whether our T-type oxy- 
sulfides are locally stable against distortion into the T'- 
type structure by performing a full structural relaxation 
in a 2x2x1 supercell with no symmetry constraints. We 
choose SC2CUO2S2 as a representative compound, since 
we believe that the smaller Sc cation will help create 
room for the large sulfur anion. The resultant structure 
(Fig. 4) shows strong octahedral rotations along the c- 
axis, denoted a°a°c~ in Glazer notation [28]. Perovskite 
structures often exhibit octahedral rotations when the 
Goldschmidt tolerance factor t — +ro)/ v2(f\B + r o) 
is less than unity as the octahedra rotate to fill the empty 
space left by substitution of small rare-earth ions. How- 
ever, we find the T-type structure to be locally stable 
and expect the rotations to be smaller in magnitude for 
the larger R =La ion. Additionally, the system remains 
tetragonal and shows no out-of-plane buckling of the Cu- 
O bonds known to be detrimental to superconductiv- 
ity [29]. We note that recent advances have allowed ex- 
perimental control over the two competing structures in 
compound synthesis [26] , which will help overcome these 
structural hurdles. 

In addition to local structural stability we check the 
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thermodynamic stability of the proposed compounds 
with respect to competing phases. We assemble the com- 
puted internal energies of formation AE = -^products — 
EVcactants in Table I (units of kJ/mol). That the quan- 
tities are positive imply the that these standard path- 
ways are not energetically favorable at zero pressure. 
However, it is known that many functional materials 
are metastable, protected from decay by large energetic 
barriers, resulting in technologically useful lifetimes [2]. 
Typical energetic deficits range from 5 to 30 kJ/mol, 
when translated to the oxysulfide compound stoichiome- 
tries. Accompanied by the fact that the volume differ- 
entials AV = ^products - Keactants are negative, high 
pressure synthesis may allow the products to form. Fur- 
thermore, we included in Table I compounds of the form 
i?2CuC>3S, where only half of the apical oxygens are re- 
placed by sulfur. The configurational entropy of the 
sulfur-oxygen arrangement in the charge-reservoir blocks 
will further encourage product formation. These chal- 
lenges in stability suggest the oxysulfides could be an 
ideal system for other materials synthesis techniques like 
molecular beam epitaxy. 

Turning to the orbital distillation theory, we show in 
Fig. 5 the effect of sulfur substitution on the energy split- 
ting AE = e x 2_ y 2 — e z 2, which quantifies the strength of 
the 3d z 2 admixture into the 3d x 2_ y 2 band. We find that 
apical sulfurs, being a larger, more covalent ion than oxy- 
gen, generally decreases AE, causing stronger mixing of 
the d z 2 orbital into the d x 2_ y 2 band, which suppresses 
Tc.max in the two-orbital model. Barring competing or- 
ders, the two proposals predict diverging effects: the re- 
duction of ed — e p should enhance T Cjmax in the oxysul- 
fides, while the increase mixing of d z 2 suppress T cmax . 

In conclusion, we have designed via first-principles a 
family of copper oxysulfides, based on the well-studied 
cuprate La2Cu04, which offer broad tunability in the 
charge-transfer energy. According to our proposal, the 
replacement of oxygen by sulfur in the apical position 
will enhance the superconducting T Ciinax , while the or- 
bital distillation proposal predicts the opposite effect. 
Although we recognize that many other mechanisms can 
drive changes in transition temperatures, we nevertheless 
believe that exploration of the oxysulfides are a valuable 
step towards testing the proposed theories. We find that 
molecular beam epitaxy is the preferred synthesis route 
since epitaxy will reduce competing orders and tip the en- 
ergetics away from competing phases. Natural extensions 
of this work include generating a detailed ab initio under- 
standing of the structural handles controlling electronic 
parameters and pursuing similar chemical engineering of 
crystal field levels in the Bi, Tl and Hg-based cuprate 
families. 
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21053 and KITP grant PHY-11-25915. G.K acknowl- 
edges support from NSF grant No. DMR-0906943. 
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FIG. 5. The energy splitting AE — E x 2_ y 2 — E z 2 between 
the onsite energies of the d x 2_ y 2 and d z 2 orbitals among the 
copper oxysulfide family (squares). For comparison is the ac- 
tual compound La2Cu04 (circle). The physical range of AE 
spanned by the single-layer cuprates is shown with a shaded 
bar (from Sakakibara, el. al.) [15]. The AE in the oxysulfides 
are generally smaller, implying stronger mixing of the d z 2 or- 
bital with the in-plane d x 2_ y 2 orbital. Thus, barring other 
competing orders, sulfur substitution should enhance T c , max 
according to the charge-transfer theory, while T c ,max will be 
suppressed according to the two-orbital theory. 
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